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Abstract 



■ ■, The linear and renormalized nonlinear kinetic theory of drift instability of plasma shear 

Qh. flow across the magnetic field, which has the Kelvin's method of shearing modes or so-called 
non-modal approach as its foundation, is developed. The developed theory proves that the 

^ ' time-dependent effect of the finite ion Larmor radius is the key effect, which is responsible for 

Q^l the suppression of drift turbulence in an inhomogeneous electric field. This effect leads to the 

^ I non-modal decrease of the frequency and growth rate of the unstable drift perturbations with 

. time. We find that turbulent scattering of the ion gyrophase is the dominant effect, which 

(— I ■ determines extremely rapid suppression of drift turbulence in shear flow. 
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I.INTRODUCTION 

^ . After seminal investigations of the regimes of the enhanced plasma confinement [1], the understanding 
• the role of flow shears in turbulent transport becomes one of the major issues of tokamak plasmas 

^ ' physics[2]. The stabilizing effect of the flow shears on various drift waves is recognized as one of the 
K*" ■ essential elements in the formation of core and edge transport barriers. A decisive step in under- 

^ ' standing the role of flow shears in the reducing of the turbulent transport and formation of transport 

5t 1 barriers has been made with the development of the gyrokinetic theory for describing turbulence in 
plasma with large shear flows [3]- [7] and the development of numerical codes[8]-[10] able to solve the 
set of gyrokinetic equations. In spite of the great progress in numerical investigations of tokamak 
turbulence, there still remain some key issues in analytical investigations of the turbulence in plasma 
shear flows, which deserve further clarification. One of the most challenging problem in the theory 
of plasma shear flows turbulence is the development of analytical methods of the investigations of 
the long-time evolution of the turbulence in shear flow governed by Vlasov-Maxwell system. Ap- 
plications of the gyrokinetic approach to analytical investigation of the instabilities in plasma shear 
flows has been made in series of papers [11]- [13], in which the electrostatic ion-temperature-gradient 
mode was considered. In all these papers, the spectral transform in time was applied and perturbed 
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electrostatic potential $ was considered in canonical modal form, $ ~ exp {—iuit). That modal form, 
however, unsatisfactory represents the long-time response of plasma shear flows across the magnetic 
field due to space-dependent Doppler shift and stretching of waves pattern by shear flows[14]. In 
fact, that approach gives results which are valid only for times limited by the condition t <^ iyo)~^ , 
where Vq is flow velocity shear. 

One of the most effective approaches to the analysis of the temporal evolution of plasma 
turbulence in shear flows is method of shearing modes or so-called non-modal approach. The essence 
of this approach, which originally was developed by Lord Kelvin[15] for fluid flows with a homogeneous 
velocity shear, consists in transforming the independent spatial variables from the laboratory frame 
to a frame convected with shear flow and studying the temporal evolution of the spatial Fourier 
modes of perturbation without any spectral expansions in time ([14], [16] and references therein). The 
transformation to the coordinates convected with shear flow eliminates the explicit spatial dependence 
related to shear flow from the convective derivative in governing fluid equations. This transformation 
not only simplifies governing equations, but is also principally indispensable. The temporal evolution 
of a separate spatial Fourier harmonic with a definite wave number can only be analyzed with 
convective coordinates; it is in contrast to the laboratory set of reference, in which spatial Fourier 
harmonics are coupled due to velocity shear. 

Kinetic effects, such as finite Larmor radius effects, Landau and cyclotron damping and the 
numerous resulting kinetic instabilities, which are naturally not involved in the fluid description of 
plasma shear flows, require the development of a non-modal kinetic description of plasma shear flows. 
Note also, that because of the shearing of perturbations in shear flow, the component of the wave 
number along the direction of the velocity shear experiences secular growth with time[14],[16], and 
therefore results obtained on the base of the fluid description are valid only for flnite times at which 
the condition k±^ {t) pi <^ 1, {k± is the component of the wave number across the magnetic fleld and 
Pi is thermal ion Larmor radius) for initially long wavelength perturbations with k± {t — to) pi <^ 1 
holds. All this requires the development of the non-modal kinetic theory, which can properly describe 
the long-time evolution of the perturbations in shear flow. Motivated by these requirements, we have 
undertaken to work out the kinetic theory of plasma shear flows, which has Kelvin's method of 
shearing modes, or the so-called non-modal approach, as its foundation. The present paper focuses 
on the development of basic equations for linear and renormalized nonlinear non- modal kinetic theory 
for electrostatic perturbations of shear flow with homogeneous velocity shear and their solutions for 
kinetic drift instability of plasma shear flow. As a flrst attempt of such analytical investigation, we 
consider here the case of the homogeneous shearless magnetic fleld. In Section II, the governing linear 
integral equation for electrostatic potential is derived. In Section III, the solution of that equation 
for drift kinetic instability, which displays the non-modal evolution with time of the electrostatic 
potential in shear flow, is obtained. The nonlinear integral equation for electrostatic potential, which 
accounted for the scattering of ions by shearing modes with random phases in plasma shear flow, is 
obtained in Section IV. The solution of that equation, which describes the suppression of kinetic drift 
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instability in shear flow is obtained in Section V. A summary of the work is given in Conclusions, 
Section IV. 



II.VLASOV-POISSON SYSTEM OF EQUATIONS 
IN SHEARED COORDINATES 

Our theory is based on the Vlasov equations, 

dF^ + + ^ f Eo (r) + i [V X B] - (r, t)) = 0, (1) 



dt dr rua \ c ' y 9v 

for velocity distribution function Fa of a species {a — i for ions and a = e for electrons) in inhomo- 
geneous electric field Eq (r), directed across the homogeneous magnetic field B, and Poisson equation 
for the perturbed electrostatic potential ip (r, t), 

A ip (r, t) = -47r ^ / (v, r, t) dVa, (2) 

a=i,e 

where is the perturbation of the equilibrium distribution function Foa, Fa — Foa + fa- The kinetic 
theory for shear flows is developed as a rule[4, 11] in a frame that is shifted by Vq (x) in velocity 
space, but is unchanged in configuration space, leaving the inhomogeneous convective terms in Vlasov 
equation. The starting point of the derivation of the basic equations of the nonmodal kinetic theory 
is the transformation of Vlasov- Poisson system to convective (sheared) coordinates in velocity and 
configuration spaces. 

We transform in Eq.(l) the variables {t, r, v) onto new spatial variables (^c, Tq, Vq), connected 
by the relations 

t^tc, V = Va + Ua (r„, tc) , r = + y U« (r^, tic) dhc, (3) 

*(0) 

or 



(4) 



Va = V - V„ (r„, t), ra = r - J \a (r, h) dh, 

Ho) 

with set of reference moving with velocity V„ (r, t) = Uq, {va, tc). In Eqs.(3),(4) t(o) denotes the time 
at which shear flow emerges. For the transformation of Eq.(l) to new variables we use the relations, 
which follows from Eqs.(3),(4), 

d_^d__ dVia(r,t) d _ _d_^ 

dt die dt dvia ' dua' 

t 

_d_^_d_ dvia^_dVia f ^Y^^-^llAdt (6) 

dvi dvia drj drj ' drj J drj ^' 

*(0) 
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With these relations Eq.(l) takes the form 

I'd? 



dtr. 



- {Via + Ui, 



t 

{ra,tc)) J 



^a: ) 

-cLti - 



Ho) 



drja 



-Via a a ^ ^ "n — 

ori ovja rrinC 



dt nia 



1 



-^(Eo^ (r) + - [Va X B], 



ea dip (r, t) \ dF {tc, 



ma dvi 



dvi, 



0, 



(7) 



where dVia/dt — dVia/dt+Vja (dVia/drj). When the velocity Via (r, t) is such as it vanishes the square 
brackets on last line of Eq.(7), Vlasov equation for shear flows, for which Uia {va, tc) dVja (r, t) /dri = 
0, contains only velocity shear parameter, V^, (instead of ^^(r) with laboratory frame variables). 



dFa ^ dFa 



dF^ 



dF„ 



dF^ 



dt 



dx 



■Jay - VaxVat) + OJ^aVay^;^ " (^^ca + V'^) '"ocx^ 



ay 



ea f dip .r,dv\ dFa dFa dip dFa dip dFa _ 

-Vj—\- \-Vaz- — - = 0. 



(8) 



nia \dx dy J dvax dza rua dy dvay rua dza dv. 

In this paper we consider stationary plasma shear flows across the magnetic field with homo- 
geneous velocity shear, for which 

c dEQ dVo 



V« = Vo 



(9) 



with Eo (r) = (dEo/dx) xe^ and dE^/dx — const. For this case the transformations (3) have a form 



t ^C) 



y ct?/ 



Van ~\~ Vr.X, 



-^ai y Ua ~l~ ^Q-^atci 



Z = Zn 



(10) 

(11) 



where Vq — dVo/dx, and t(o) =0 was assumed. 
With leading center coordinates. 



■^a -^a 



sin (01 - ^/Jl^UJcat) , 



^a Va 



COS (01 - y/]j^UJcat) - V^t {Xa - Xa) , Zi = Z - V^t, 



and velocity space coordinates (see, also Ref.[19]), 

Vax = V± COS 0, Vay = Siu 0, = 0i - ./Jl^UJcat-, = Vo 



(12) 



(13) 



where /Iq, = 1 + Vq/Uco > 0, the equation for the perturbation fa of the equilibrium distribution Fqq, 
function {Fa = Foa + fa) has a simple form. 



+ 



dfa 

dt ' ma^/lI^UJca \dXadYa dYadXa 

dip dF^ 



dip dfa dip dfa \ _^ eg ^fiQjJca ( dip dfa dip dfa 



nia 



i)i dv^ 



Ca dp dfa _ ea 



rUa dza dv 

za 



1 



Oa 



dYa dXa 
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^(jJca dp dFoa _^ dp dFoa 



V± d(pa dv±a dZa dv^ 



(14) 



We assume that jV^'/'^cal ^ 1 and put in what follows /la — It is interesting to note, that the 
equilibrium distribution function Fq^, which in laboratory frame contains the spatial inhomogeneity 
resulted from electric field Eq (r), does not contain such inhomogeneity in convective coordinates (see 
Appendix 1). In what follows we consider the equilibrium distribution function Fjo as a Maxwellian, 

exp f-'lpi] . (15) 



Foa — 



assuming the inhomogeneity of the density of plasma shear flow species on coordinate Xq,. It follows 
from Eq.(8), as well as from Eq.(14), that for Vq = const, these equations do not contain the spatial 
inhomogeneities, originated from inhomogeneity of the flow velocity Vq (x). Therefore, the spatially 
homogeneous, but time dependent, Eqs.(8) and (14) may be Fourier transformed over the variables 
Xa,ya,Za with Conjugate wave numbers k^, ky and kz and the temporal evolution of the separate 
spatial Fourier mode of the perturbations of the distribution function, fa and of the electrostatic 
potential, ip (k, t) may be traced upon. On this way we present the potential ip (r, t) in form 

{^a, Ua, Za, t) — J" </? {kx, ky, kz, t) 6 ^ ""^ y2/a+ z °' (JJ^^dl^yd]^^ 



j {kx, ky, kz, t) exp [ikxX^ + ikyY^ + ikzZ^ 
k± {t) v± 



—I 



sin (0 - ujcat -0{t)) 



where 



kl {t) = {kx - V^tkyf + kl 



(16) 



(17) 



and tan^ = ky/{kx — V^tky). The solution to linearized Eq.(14) is calculated easily for any values of 
the velocity shear rate Vq and it is equal to 

t 

oo oo „ 

fa{t,kx,ky,kz,v±,(f),Vz,Zi) = — ^ dtiip{ti,kx,ky,kz) 



n=—oo ni=— oo 



X exp ^ — ikzVz {t — ti) + in (0i — coct — 9 {t)) — ini (0i — cocti — 9 {t)) ^ 



/C_L (t) 



J, 



ni 



kl. {ti) v± 



ky dFa , UcnidFa , , dF^ 
H \-ki 



UcadXa V± dvj 



^ dvz 



^fa {t — to, kx, ky, kz, V±_4), Vz) ■ 



(18) 



The Poisson equation for separate spatial Fourier harmonic ip (k, t) gives the governing integral 
equation [17] 

t 

oo „ 

{kx - V^tky)' + kl + kl] ip{k,t)=j2^ J2 / ^^l'^ 



to 



X 



In {k± {t) k± {ti) pI) e"^''"(''iW+^i(*i))e"5'''"Tc.(*-*i)''-»""<=«(t-ti)-m('?W-e(ti)) 
X [kyVda - nuca + iklv^a {t " ^i)] - 47r ^ eaSna (k, t, to) , 



(19) 



where 



47r Ca^Ua (k, t, to) 



a=i,e 
oo 



rv=.i. p _ ^ ' 



(20) 



and where it = to, k, v±, v^) is the initial, determined at i = to perturbation, assumed here as not 
dependent on 0, of the distribution function F^. It follows from Eq.(19), that initial perturbation 
(p(k,t = to) of the self-consistent electrostatic potential is equal to 



(/? (k, t = to) 



kl {to) + kl 



ea^Ua {k,to,to) . 



(21) 



Eq.(19) presents minimal model of the linear kinetic theory of plasma shear flows, which incorporates 
the shear flow effects. In Eq.(19) these effects are concentrated in the time dependent arguments of 
the Bessel function and manifest itself as a time-dependent finite- Larmor-radius effects. 

The integration by parts of the first term on the right of Eq.(19) gives the integral equation, 
which appears to be more convenient and transparent for further analysis. 



{k. - V^tkyj' + kl + kl ip{k,t) = -J2 A^'^ ^) + E E / ^^^^ i^fi ^1 



dti 



to 



x/„ {k^ it) k^ {ti) pI) e-i''«('=iW+'=i(*0)— e^W-'^C*!))! e-l'^.^- 



A p ' 



to 



X exp 



-IpI {kl (t) + kl (ii)) - in {9 (t) - e (t,)) - hlvl^ {t - t^f - inw,^ {t - t,) 
-47r ejria (k, t, to) + ^ -^(p (k, to) Pa {t, to) , 



a=i,e 



(22) 



where v^a — cT^/eBLn is the diamagnetic drift velocity, L^^ — —d\n.noa{x) jdx, Pa is thermal 
Larmor radius, and 



Pa {t, to) = E (^^ (^o) ) 



X exp 



—pI {kl (t) + kl (to)) - in {0 (t) - (to)) - ^ky^^ {t - to)'' - inouca {t - to) 



(23) 



Note, that Pa {to, to) = 1. 

It is important to note the alternative, explicitly causal representation of Eq.(22) with function 
$ (k, t) — (fi (k, t)Q{t — to), where Q {t — to) is the unit-step Heaviside function (it is equal to zero 
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for t <to and equal to unity for t >to). That equation for $ (k, t) has a form 

{K - V^tkyf + kl + kl $ (k, t) + 47r ^ eju^ (k, t, to) 

a=i,e 

t ^ t 

a=i,e a=i,e n=— oo ^ 

oo t 



to 



X exp 



■^pI {kl it) + /cl (ti)) - (t) - 9 ih)) - ^ky^^ {t - Uf - iTiu:,^ [t - h) 



(24) 



Using the quasi neutrahty approximation with {k^ (t) + k'^) A|,^ <^ 1, and averaging this equation 
over the time t 3> u~^, we obtain from Eq.(24) the equation, which is relevant for the analysis of the 
low frequency drift type perturbations, 



*(k,ti) 



1 + r) - /o [k^ (t) k^ (t,) p^) e-^^K^iW+'^K^O)] I 
-i J dh^ik,t,)kyV,Jo{k± {t)k^ it,)p^) e-^''?K(*H'=i(*i)) 

to 

t 

-Jdh-^ ($ (k, t,) lo {k^ (t) k^ (h) p^) e-^^?('=iW+'=i(*0)) _ 

to 

t 

i J dh<^> (k, h) kyVaJo {k± it) ki_ ih) pI) e-i''?('=i(*)+'=i(*i)) (^^ _ ^-^^klvUt-ti?^ 



to 



+r / dti 



(25) 



where r = Tj/Tg. Eq.(25) introduces obvious time scales which determine different stages of the 
temporal evolution of the electrostatic potential in plasma shear flows. As it follows from Eq.(17), 
the time dependence of k^_{t) is negligible in times t <C (V^')~^ ^^"^ ^e neglected; on that 
stage Eq.(25) determines the ordinary modal evolution of perturbations as in plasma without shear 
flow. In times (VJ)"^ <^ t <^ = {yikyPi)~^ time dependence of k±_ gradually enhances the non- 
modal evolution of the potential with time. In times t tg, non-modality becomes the dominant 
effect of the temporal evolution of the potential $ (k, t). In next Section we will obtain by successive 

approximations the approximate solution to Eq.(25) for long wavelength, k± (to) Pi < 1, perturbations 

1 _ e-lfe^li(t-to)' 



with weak ion Landau damping, for which 



<^ 1. The solution will be obtained 
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without application the spectral transformation over time for the case of a weak velocity shear, or a 
small time, for which condition |V^'| < t <^ ts is met, and for long times, t ^ tg- 

III.LINEAR NON-MODAL EVOLUTION OF THE 
KINETIC DRIFT INSTABILITY OF 
PLASMA SHEAR FLOW 

If k± (t) Pi < 1 at time t — i(o) = at which the shear flow emerge , i.e. Vq = VqQ (t), we will get 
k± {t) Pi < 1 on times t < tg throughout. By using the approximation 

+ (^KkyV^pj {t + 1,) - hy, (vj)' {t' + 1?)) (^) (26) 

in Eq.(25), we present Eq.(25) in the form 



to 

t 



t 

f , f (k, ti) . „ , ^ 



to 





t 



ky (ijV^'t^ 2ajt^ 






+ ^ / [^^^ + ^^.-.e$ (k, to) (l - e-^'=J^-(-*o)^) , (27) 

to 

where bi — 1 — k\pl, = r + k\p1 and 

a; (k) = -- Vdi- (28) 

The first term in the right side of Eq.(27) determines the ion Landau damping; this term is the same 
as in plasma without shear flow. The next three terms originated from shear flow and determine 
the corrections to the frequency and ion Landau damping provided by shear flow. The right side of 
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Eq.(27) is small for (V) ^ < t < tg, r < 1 and for weak ion Landau damping. Therefore the solution 
to Eq.(27) we seek in the form 



$ (k, t) = C exp {-iu (k) t + iu (k, t)) . 



(29) 



Inserting Eq.(29) into Eq.(27) and neglecting the derivative dp (k, t) /dt in the right side of Eq.(27), 
we obtain for u (k, t) following equation, assuming that to — ^ — oo, 



dti^ {k,ti 



to- 



dh 



iu (k) ti 1 + 



a,; 



-U 



0, 



(30) 



where 



.uj (k) (g; (k) - kyV^i) k 



2 \^/2kzVT^ 



^ .^ U (k) (U (k) - kyVde) [^^rf ^i^) 



(31) 



with uj (k) defined by Eq.(28) and W {z) — exp(— 2;^) ^1 + J dC,e'^ j. It is interesting to note, 

that Imt^o; (k) = 7 (k) is right well known growth rate of the kinetic drift instability [18] determined 
in plasma theory approximately as 7 (k) = — Im£ (k, u; (k)) / (9Re£ (k, u (k)) /du (k)), where e (k, u) 
is ordinary electrostatic dielectric permittivity of the plasma in magnetic field. By equation to zero 
the expression in square brackets in Eq.(30), we obtain simple solution for v (k, i), which determines 
the non-modal evolution of the potential $ with time. 



$ (k, t) = $0 exp 



-iu: (k) i ( 1 - ^7^1 + iRe5u (k) i + ("7 (k) - ^ 



2aitl 



(32) 



As it follows from Eq.(32), non-modal effects, which reveal in non-modal reduction of the frequency 
and growth rate, are negligible aXt <^ts and become dominant at i ~ ig. Note, that for r » k\pl the 
time (I'/'^ts is approximately equal to time ^2 = iyok^Ps)'^ of the transition to strongly non- modal 
regime in the fluid theory of the drift turbulence of the plasma shear flow[16]. 

In the laboratory frame of reference spatial Fourier mode (32) is observed as sheared mode 
with time dependent component of the wave number kxiiah) — kx — kyVgt directed along the velocity 
shear, and therefore is quite different from the normal mode assumption. 



* (r, t) =jik.j ik. j exp + ^k.y + - V.'t.) 



X exp 



-icu (k) t ( 1 - ^-^^] + iRedcu (k) t + (^7 (k) - ^ 1 1 



ttihi 3i? 



(33) 



The mode shearing, which is other non-modal effect, becomes pronounced in times t > 7 ^ when 
Vq ~ 7. Only in times, for which |V^'t| ^ 1 solution (33) has a normal mode form, $(k, t) ~ 
exp yik^x ~\~ xkyU -\- ik^ z — iu) (k) t). Because of the time dependence kx(iab) = kx — kyV^t of the wave 



9 



number component along the flow shear, the modes in the laboratory frame in times t ^ {Vo)~^ 
become increasingly one-dimensional zonal-like as the perturbed Ex B velocity tilts more and more 
closely parallel to y-axis. 

The exceptional advantage of the non-modal approach is a possibility to perform the analysis 
of a plasma evolution on any finite time domain. Now we consider the temporal evolution of drift 
perturbations at times t > Iq ^ tg on the base of Eq.(22) averaged over the time t ^ u~^. Note, 
that for the averaged Eq.(22) , 



Pi (t, to) = lo {k± (t) (to) pI) exp 



-\pnkl{t) + kl{to))-\klvUt-tof 



(34) 



and Pi (fo, to) — lo (^1 (io) p}) e~''<*^J-(*°). We consider the times t > to 3> ts, for which k±_ (t) becomes 
large enough, so that k± (t) pi > k± (to) Pi > 1 and Pi (to, to) ^ (V^k± (to) pi) ^. The initial value 
(p(\s.,t — to) for the averaged potential in that case have to be corrected and is equal to 



(k, i = to) = (p (k, to) 



27rto 



where ip (k, to) is the initial value for Eq.(22) without averaging. For the times considered 



lo {k± (t) k^ (ti) p2) e-i^?('=i(*)+'=i(*0)-^fe?4.(*-*i) 



y/27rtri 



exp (t - tiY 



where 



«i = 7^ + k,n 



With this approximation, Eq.(22) becomes 



(35) 



(36) 



(37) 



t 

(1+.)/ 



dti 



dtp (k, ti 
dti 



dti 



d 

dti 



to 



to 



(^(k, ti) 



(27rtti)'/' 



e 



+ikyVdi / dtiip(k,ti) 



to 



(2'Ktti) 



1/2 



e 2 



f(t-ti)' 



+r| dti (^^^ + .V.e^(k,t0) e-i'=?^-(*-*^)= 

to 



+V (k, to) 



(2Tltto) 



1/2 



e 2 



tit-tof 



to 



1/2' 



(38) 



In Eq.(38), electron term may be omitted, because for times t — to > tg and V^' ~ 7 = kyVdik^p'^, 
where 7 is the growth rate of the drift kinetic instability. 



1 



exp ( (t - tif ) < exp ( —klvlf 



exp 



kl rrij k\Ll 
Tkl me (k^pi)'' 



< 1. 



(39) 
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t 

In zero approximation over ts/t we have for ip(k.,t) equation f dtid(p (kjti) /dti — with solution 

to 

(f (k, ti) /dti = ifo = const. Accounting for in first approximation small right hand side of Eq.(38), 
we seek solution for 93 (k, t) in the form 

(/?(k,i) = (^oexp(i/(k,i)), (40) 

where v i\s-,t) — 0{ts/t). For times t > to ^ k,"^ we can omit in Eq.(38) the terms with small 
exponents exp (— (fi;^/2) {t — to)^) and obtain the following solution for cp (k, t): 

which gradually becomes a zero-frequency cell-like perturbation. 



IV.NONLINEAR PHASE SHIFT IN SHEAR FLOW 

Nonlinear terms in the left hand side of Eq.(14) can drastically change linear non- modal solutions 
(33) and (41). Because of the nonlincaritics, accounted in the left side of Eq.(14), variables X, Y, v±, <p 

and z become coupled and determined by equations of characteristics, 

dX dY dv± d(j)i dv^ 



dt = 



dip e dip e Ud dip e cUd dip e dip 



rriiUJci dY rriiUJci dX rrii w_l dcpi rrii v± dv± rrii dzi 

^ dfi 

e dip dFio e ujd dip dFio e dip dFio ' 



(42) 



rriiUJci dY dX rrii ^ ± d(j)i dv± rrii dz-y dv^ 
Last equation in system (42) gives nonlinear solution for the perturbation of the ion distribution 
function /j with known Fjo, 

t 

f.^ ^ [ "J_^^^ _ cUddi^dFio _^ dip_dFia 
' m J cuci dY dX v± d(/)i dv± dzi dvz 

in which coordinates X, Y, v±, are X — X + 5X, Y = Y + 6Y, (j) = (j) + 6(f), where X, Y, are the 

guiding center coordinates averaged over the turbulent pulsations, and SX{t), SY{t), 5(f) are random 

ion orbit disturbances due to their scattering by electrostatic low frequency drift turbulence. The 

disturbances arc assumed sufficiently small and, after the averaging over the times t S> (cUci)^^, they 

are determined by the equations 

t t 

5X^-^1 %dt, = [ dt, [ dk^ (k, to k,Jo e^*, (44) 

rriiUJci J dY B J J \ J 

to to 

SY = ^'f = -I / rfii / dk^ (k, (,) Jo ( '^^^) e», (45) 

rriiUJci J dX B J J \ cjd J 

to to 

t t 



= ! = ^fdtJ dk„ (k,tO (*,) 7, ( ^i<M£i) e-\ (46) 

rrii v± J dv± mv± J J \ J 



5(f) 

rrii vi. 

to to 
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and Sv± = 0. In Eqs.(44)-(46), * = k^X + kyY + k^z + k{ti) Sr {ti), and ik {t) Sr {t) denotes the 
phase shift resulted from perturbations of the ions orbits due to random ion-waves interactions, 



k (t) 6r (t) = kJX (t) + ky6Y (t) - ^^ ^^^ cos (0 - 6) 50 (t) 



(47) 



In Eq.(47), the scattering of ions along the magnetic field is ignored. In variables X, Y, 0, v±, 
the averaged over the time t ^ to^^ and over initial phases of the drift perturbations solution for 
fa {t, ky, kzi vx, 0, Zi) in drift frequency range has a form 

t 

fa(t,kj;,ky,kz,v±,(l),V;„zi) ^i— / dti(p{ti,kx,ky,kz) 
X exp ( - tk,v, {t -h)-^ ((k (t) Sr (t) - k (ti) Sr {t,)f) ) 



UJc. 



ky dFg ^ dF^ 

UJca dXa dVz 



+fa {t = to, kx, ky, kz, V±(p, Vz) , 



(48) 



in which fa {t — to, k, v±, v^) is the initial, determined at f = to perturbation, assumed as independent 
on (f), of the distribution function Fa- In Eq.(48) Gaussian distribution for ions orbit disturbances is 
assumed. We use Eq.(48) in Poisson equation for the potential (p (r^t). 



A (/? (r, t) = -47r ^ Ca J fa (v, r, t) dVa, 



(49) 



and obtain integral equation for separate spatial Fourier harmonic ip (k, t) for electrostatic potential, 
in which effect of the turbulent scattering of ions on sheared drift modes of random phases is accounted 
for, 

t oo 

k'^ {t)ip{k,t) = ^ ^2 ^,,2 / dti(p{k,ti) J dv±v±exp 



^Da^Ta 



to 



'Ta 



X 



) exp (-^^Xa it - tif -\{{k it) Sr it) - k it,) Sr (iO)'>) 
[kyVda + iklv^a {t - ti)] ~ 47r ^ CaSUa (k, t, to) . (50) 



It follows from Eqs.(44)-(46), that for k^ ~ ky we have the estimates SX ~ SY and 

kyJo 



k JX (t) 



it) —S(l> (t) 



^2 7 ik±{ti)vx_ 



kl (ti) Ji 



ki. (ti) w± 



(51) 



At initial times of the evolution, t <^ (V^') ^ , for long wavelength perturbations with k± (t) v± <^ Ud, 
the nonlinear phase shift is determined mainly by SX and SY. For short wavelength perturbations 
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with k± {t) v± > LOci for these times the if) term is important as well. At these times non-modal 
effects are negligible. At times t > {Vq)~^ for k± (t) v± <^ ood 



kJX (t) 



and for k± {t) v± > oud, 



k± (t) —5<l){t) 



kJX {t) 



1 1 



3 ' 



(52) 



k± {t)—S<P{t) 



2 • 



(53) 



It follows from Eqs.(52), (53), that at times t > {Vq)~^ turbulent scattering of the angle (50 (t) is 
the dominant process in the formation the turbulent shift of the phase of the electrostatic potential. 
Now consider the average phase shift term for times (V^kyPi)'^ > t > {Vq)~^, but for k±v±/oJci < 1. 
In that case k± (t) f=:i kyV^t and 



Ji 



k± (t) v± \ k± (t) v± kyVQtv± 



IjJr 



2Ur 



IjJr 



(54) 



In convective coordinates, solution for potential ip (k, t) at kyV^tpi < 1 (and V^t > 1) has a modal 
form 



(k,i) = <^(k,to)e^"^'^)*+^('^)*. 



(55) 



where u (k) and 7 (k) are the frequency and growth rate of the kinetic drift instability. For times 



t > (Vq) ^, the approximation 

{{k{t)5r{t)-k{t^)Sr{t,)f) 



ct 



2ul 



t t 



j dt[ J dt'^ j dki |(^ (ki, to) I' kiy exp [7 (ki) (i; + Q + iu (ki) {t\ - Q ] 



X {t[t'^fexp 



~{k{t[){6r{t[)-Sr{t',)f 



(56) 



is valid. With time variables r — ti — t2 and t — {ti + ^2) /2, Eq.(56) becomes 



2a; 



((k(i) {6cl>{t)-5cf>{t,))f) 



J dr J di+ J dr j di 
\i-t ti-5 ti+i J 



13 



X 



J dk,\ip (ki, to) I' A;t^t4g27(ki)t+ia.(ki)r 



/ t+|(T+ri) 



X 



t+i(T-Tl) \ 



dri J dti + J dri J dti 

\-r f-i(r+n) f-l(r-ri) / 

X exp (27 (ks) ii + iuj (ks) Ti 

■i ( (k k + in) (& (t; + ir.) - Sv U - \n 



d\^2 \^ {\^2,tQ)f k2yi{ 



(57) 



The integration over ii is performed over narrow interval (t + | (r ± ri) , t — | (r ± ri)) . In such a 
case the integrals over ti may be approximately calculated as 



t+i(r±ri) 



(58) 



f-i(r±ri) 



and for Eq.(57) we obtain 



2ui 



(k(i) (50(t)-50(ii)))'> 



*+§ i-ii 
dr j di+ [ dr [ d.t 

ti-t ti-r 



ti+i 



exp I — T 



y" dki |(/p(ki,to)|' A;yV^('^i)*+^'^('^^)^ 

cik2A;2\|^(k2,to)|'e2^('^^)*: 



a;Mk2) 7' 

where C (k2, t) resulted from infinite sequences of inner integration in exponential of the 
k (ti + ^Ti) (5r (ti + |rj) — 5r (ti — over ii and determined by integral equation 

^ ^^'^'^ ^ Ib'J I ^^^^^^ l^(k2,^o)r e^-C^^)*-^^. 



(59) 



a;2 k. 



(60) 



Changing r — )■ — r in first integral over r in Eq.(59), we perform the integration over r and obtain 
simple result 



2a;2. 

Cl 



{{k{t){S(j){t)-6(j){ti))f) = 2t^ / di 



;C{k,i) 



(61) 



V.NONLINEAR NON-MODAL EVOLUTION 
OF THE KINETIC DRIFT INSTABILITY 
IN SHEAR FLOW 

Now we use Eq.(61) with v\ changed on v'^^ in Eq.(50) and integrate (50) over v±. We obtain in 

t 

ion terms additional multiplier, exp(— J dtt~^C (k, i)), as compared with linear non-renormalised 

ti 
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equation (22), 



e (t) (k, t) 



to 



X 



X 



[/cj^Vde + iklv^^ {t - ti)] - 47r ^ CaSUa (k, t, to) • (62) 



For the perturbations of drift type, the equation for $ (k, t) = (p (k, t) 6 (k, t) has a form 

t 

J dt, { A$ (k, t,) (1 + r) - [A (k, t,) /o (A;^ (t) k^ ih) pI) e-^^f 



to 



-J- ( $(k,tO ( 1 /o(fcx(t)fcx(ti)p.^)e-^''K'=i(*H^i(*0) 



; j dh^ (k, ti) A;,^;,,7o {k± (t) k^ (h) pi) e-^''K'=iW+'=i(*0) 



to 



-i y (k, t,) kyvajo {ki. it) k^ ih) pI) e-i''?('=iW+'=i(*0) | i _ e 4 ^ 



+ 



to 



to 



(63) 



If kj^pi < 1 at time t = (at which the shear flow emerge), we will get k±_ {t) pi < 1 on times t <ts 
throughout. By using the approximation 



k + (^k^kyV^ {t + h) - hi {v;f {t' + tl)^ pIq it) , 



(64) 



where bi ~ 1 — k^pf, and 6 {t) indicates that the shear flow emerge at t = 0, we present Eq.(63) in 
the form 

t 



/*( 

to 



d$ (k,ti) 



+ iuj (k)<l>(k,ti) 



= / dti 

a. 



j + (k, io) (l - e-f ^*('-'.'') 

to 



/ (^^^ + ikyva^ (k, ti)) [ 1 - exp j 

to \ L ti 



C (k, t) 



t2 
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+^ / *i [^^^ + (k. (l - , (65) 

to 

where u; (k) is determined by Eq.(28). The right hand side of Eq.(65) is small for {Vq)~^ < t < tg, 
for T < 1 and for weak ion Landau damping. Therefore, the solution to Eq.(65) we seek in the form 

$ (k, t) = C exp {-icu (k) t + iu (k, t)) . (66) 

Applying the procedure of the solution of the integral equation (25) to the renormalized version of 
that equation, (64), we obtain for {Vq)~^ < t < tg the solution to Eq.(65) in the form 

iu (k) t (l - t + iReSu (k) t 



$ (k, t) = $0 exp 



(67) 



where C (k, t) is determined by the equation 

CiKt)-^,klp^^^ I cik,|^(k„t)|^C(k„t)-^. (68) 

If we omit linear non-modal terms in Eq.(67), the condition of the balance of the linear modal growth 
of the kinetic drift instability and non-linear non-modal dumping is determined by the equation 
7 (k) = C (k, t). By using this equation in Eq.(68), we obtain the equation, which determines the 
time, at which that balance occurs, 

^-8i^SP?/*l.(>'..')lS(k.)^. (69) 

The effect of the shear flow reveals in the reducing with time as {V^ty^ the magnitude of the 
growth rate in the left part of the balance equation (69). That causes rapid suppression of the drift 
turbulence. The evolution of drift turbulence in times t >ts continues as strongly non-modal process, 
for which Markovian approximation, which is admissible for the solution Eq.(23) with small growth 
rate and non-modal terms with respect to the frequency cu (k), ceases to be valid. 



VI.CONCLUSIONS 

In this paper, by using method of shearing modes or non-modal approach, originally developed by 
Lord Kelvin[15] for fluid descriptions of fluid shear flows, we develop for the flrst time non-modal 
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kinetic theory of plasma shear flow directed across the magnetic field. We obtain linear, (22), and 
renormalized nonlinear, (62), governing integral equations for the perturbed electrostatic potential. 
By using these equations we obtain linear, (32), (41), and renormalized nonlinear, (67), initial value 
problems solutions for kinetic drift instability of plasma shear flow. Obtained solutions display two 
distinct non-modal effects, which are observed in laboratory frame of reference. 

The flrst effect displays the inhomogeneous Doppler shift, which is presented by the term 
V^tkyX in exponential of Eq.(33). This term displays the shearing of waves patterns in shear flow, 
observed in laboratory frame; in time t > {kyVQ/kx)~^, turbulence becomes almost one-dimensional 
in plane across the magnetic fleld and directed almost along the shear flow. It is obvious that this 
Doppler shift is irrelevant to the suppression of drift turbulence by shear flow. 

Second non-modal effect is of principal importance for turbulence evolution in plasma shear 
flows. It reveals as a time dependent flnite Larmor radius term in governed equations (22) and 
(62). The time dependence originates from the dot product of the time dependent coordinates of ion 
gyration in sheared coordinates, (Eqs.(12)), and wave number, which is time independent in these 
coordinates. That term completely conserves its form after the transformation to the laboratory 
frame variables, in which ion Larmor orbit is almost circular in shear flow with IV^qI ^ Ud [19]. In 
laboratory frame of references, this non-modal effect is seemed as resulted from the coupling of the ion 
gyration and temporal variation of the wave number, k±(iab) = {ky + {kx — kyVQt)^) ^ of the shearing 
mode. Linear theory reveals (sec Eqs.(32),(41)) that on the times (Vq)~^ < t < (Vf^kyPi)'^ this time 
dependence resulted in the reducing the frequency and growth rate of the drift kinetic instability 
and to gradual suppression of the instability. It is important to note, that this effect is absent in the 
theory grounded on drift kinetic equation, in which effects of the finite Larmor radius are ignored. 
It is interesting to note that similar effect of the non-modal evolution of the perturbed electrostatic 
potential, which consists in reducing with time the frequency and growth rate, was discovered[14] 
in the investigations of the resistive drift instabihty on the base of the Hasegawa-Wakatani system. 
That effect for resistive drift instability originates from the time dependent polarization drift which 
is in fact the manifestation of the time dependent effect of the finite ion Larmor radius. 

Decisive impact on the temporal evolution of the kinetic drift instability has non-linear non- 
modal effect of the turbulent scattering of ions by the ensemble of sheared waves. We find that 
turbulent scattering of the gyrophase of ion Larmor orbit is the dominant effect, which determines 
extremely rapid suppression of drift turbulence by flow shear. 
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Appendix 1. Equilibrium distribution function Fq^- 



The equilibrium distribution function Fq^ {x,Vx,Vy,Vz) is governed by the equation 

Vx-K \- I Eq [x) + UcaVy I ^ OJcaVx^ = 0. [A.l) 



Prom characteristic equations 

dx dvr dv. 



y 



g 

'^x — Eq ix) + (jJcaVy '^coVx 

m 



{A.2) 



the relation 



g 

UIcEq {x) dx + UcaVydVy + UcaVxdVx — d (uJcH) — (^-3) 

follows in which is Hamiltonian of a particle in electric and magnetic fields. By using the expansions 

Eo (x) = Eo (X^) + K {Xa) {x - Xa) , (^-4) 

and 

Vy = Vy^ + Vo (XJ + {x - X^) , (^.5) 

and accounting for that Vq {Xa) — —cEq {X^) /Bq and dvy — —oUcadx, we obtain that 

-Uol^VaydXa + OJca^xadv^a = d {uJcH) = 0, (^.6) 

where identities x — Xa and = v^a were used. It follows from (^.5) that dVy = dvy^ + {Xa) dx^- 
Therefore dx^ — dx — —dvy/ujca = ~ {dvya + Vq (Xa) dxa) /(^ca and dxa — —dvya/ iJ^uJca- With 
Eqs.(13) it follows that 

^d{ujcavl) =d{ujctl) = 0. (A7) 

The same conclusion about absence of the spatial dependence in Hamiltonian in sheared coordinates 
follows from the analysis of the characteristics of Eq.(8). 

References 

[1] F. Wagner, G. Becker, K. Behringer, D. Campbell, A. Eberhagen, W. Engelhardt, G. Fussmann, 
O. Gehre, J. Gernhardt, G.V. Gierke, G. Haas, M. Huang, F. Karger, M. Keilhacker, O. Kluebner, 
M. Kornherr, K. Lackner, G. Lisitano, G.G. Lister, H.M. Mayer, D. Meisel, E.R. Mueller, H. 
Murmann, H. Niedermeyer, W. Poschenrieder, H. Rapp, H. Rohr, F. Schneider, G. Siller, E. 
Speth, A. Staebler, K.H. Steuer, G. Venus, O. VoUmer, and Z. Yue, Phys. Rev. Lett. 49, 1408 
(1982) 

[2] K.H. Burrell, Phys. Plasmas 4, 1499(1997) 
[3] A.J. Brizard, Phys. Plasmas 2, 459, (1995) 



18 



[4] T.S. Hahm, Phys. Plasmas 3, 4658, (1996) 

[5] H. Qin, Contrib. Plasma Phys. 46, 477, (2006) 

[6] G. Kawamura, and A. Fukuyama, Phys. Plasmas 15, 042304 (2008) 

[7] Lu Wang, T.S. Hahm, Phys. Plasmas 17, 082304 (2010) 

[8] R.E.Waltz, G.D.Kerbel, and J. Milovich, Phys. Plasmas 1, 2229 (1994). 

[9] R.E. Waltz, G.D. Kerbel, J. Milovich, and G.W.Hammett, Phys. Plasmas 2, 2408 (1995). 

[10] R.E. Waltz, R.L. Dewar, and X. Garbet, Phys. Plasmas 5, 1784 (1998) 

[11] M. Artun, W.M. Tang, Phys. Fluids B4, 1102, (1992); M. Artun, J.V.W. Reynders, W.M. Tang, 
Phys. Fluids B5, 4072, (1992) 

[12] G. Rewoldt, M.A. Beer, M.S. Chance, T.S.Hahm, Z.Lin, and W.M.Tang, Phys. Plasmas 5, 1815 
(1998) 

[13] M. Maccio, J. Vaclavik, and L. Villard, Phys. Plasmas 8, 895 (2001) 

[14] V.S. Mikhailenko, V.V.Mikhailenko, K.N. Stepanov, Phys. Plasmas 7, 94 (2000) 

[15] Lord Kelvin (W. Thomson), Philos. Mag. 24, 188 (1887) 

[16] V.S. Mikhailenko, V.V. Mikhailenko, K.N. Stepanov, Plasma Phys ControU Fusion 52, 055007, 
(2010) 

[17] V.S. Mikhailenko, V.V. Mikhailenko, K.N. Stepanov, Plasma Fusion Res. 5, S2015, (2010) 

[18] B.B. Kadomtsev, Plasma Turbulence, Academic, New York, (1965) 

[19] K.C. Shaing, A.Y. Aydemir, R.D. Hazeltine, Phys. Plasmas 5, 3680, (1998) 



19 



